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Electromagnetic dipole absorption cross-sections of transitional nuclei with large-amplitude shape 
fluctuations are calculated in a microscopic way by introducing the concept of Instantaneous Shape 
Sampling. The concept bases on the slow shape dynamics as compared to the fast dipole vibrations. 
The elctromagnetic dipole strength is calculated by means of RPA for the instantaneous shapes, 
the probability of which is obtained by means of IBA. Very good agreement with the experimental 
absorption cross sections near the nucleon emission threshold is obtained. 
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Photo-nuclear processes, as the absorption of a photon 
inducing the emission of a neutron or the emission of a 
photon after neutron absorption, are key elements in vari- 
ous astrophysical scenarios, like supernovae explosions or 
7-ray bursts, as well as in simulations for nuclear tech- 
nology. For a quantitative description of the relevant nu- 
clear reactions one needs to know the y-absorption cross 
section and the reemission probability, being determined 
by the dipole strength function. Direct measurements of 
the strength function in the relevant energy range (typi- 
cally 6-10 McV in medium-heavy nuclei) are not possible 
nowadays for most of the unstable nuclei passed in violent 
stellar events. Theoretical models that provide reliable 
prediction of the dipole strength function are therefore 
of utter importance. Aside from the astrophysical ap- 
plications understanding the mechanisms that determine 
the structure of the dipole strength function in this en- 
ergy region is a challenge of its own to nuclear theory. 
The present Letter proposes and tests a new approach, 
which we call Instantaneous Shape Sampling (ISS). It 
combines the microscopic Random Phase Approximation 
(RPA) for dipole excitations with the phenomenological 
Interacting Boson Approximation (IBA) for a dynamical 
treatment of the nuclear shape. ISS allows one to calcu- 
late the dipole strength function of the many transitional 
nuclei ranging between the regions of spherical and well 
deformed shape. 

Traditionally one employs phenomenological expres- 
sions for the dipole strength function [l| , which are based 
on the classical model of a damped collective Giant 
Dipole Resonance (GDR). The photo-absorption cross 
section ct 7 of the GDR is approximated by a Lorentzian 
curve 0, 0, [3, El, which may include corrections for 
nuclear deformation The damping width of the 

Lorentzian is treated as a parameter that is adjusted to 
the experiment. However, the available data have not yet 
allowed stringent tests of this extrapolation toward the 
low-energy tail of the GDR [4j. Therefore microscopic 
approaches that treat at least a substantial part of the 
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FIG. 1: (Color online) Probability distributions of the in- 
stantaneous nuclear shapes over the /3 — 7 plane. Families of 
coexisting shapes are distinguished by their color. 

damping explicitly promise improved predictive power. 

RPA [6] is the standard microscopic approach to the 
dipole strength function. It takes the Landau fragmen- 
tation into account, which describes the coupling of the 
collective dipole vibration to the particle-hole excitations. 
It also describes the splitting of the GDR caused by a 
static deformation of the mean field. However both ef- 
fects account only for a fraction of the observed width of 
the GDR. In order to obtain the experimental damping 
width of the GDR one must include, at least, the cou- 
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TABLE I: IB A parameters f, \, e_g, and equilibrium deforma- 
tion parameters f3, 7 calculated by means of the micro-macro 
method. In case of shape coexistence, two sets are listed. 
Their respective proportion in the ground state is given in 
percentage. 
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pling to the two-particle two-hole states [6J]. One may 
distinguish between two types of such states: the com- 
bination of the GDR with incoherent particle-hole exci- 
tations and its combination with coherent collective ex- 
citations. The coupling to the former is analogue to the 
collisional damping of Fermi liquids. The latter has been 
studied for spherical nuclei by means of particle-phonon 
coupling models, such as QPM 0, H, QTBA i, and 
QRPA-PC [10], which, however, meet principle problems 
in transitional nuclei. 

We suggest an alternative approach. Out of the huge 
space of collective excitations coupling to the GDR we 
select only the low-energy collective quadrupole excita- 
tions. They represent the softest mode, which couples 
most strongly to the dipole mode. The coupling, coher- 
ently and incoherently, to the other excitations is taken 
into account by averaging with a Lorentzian of the type 
that arises from collisional damping 0, 0]. In the fol- 
lowing we simply refer to this as "collisional damping" 
(CD). The typical energies of the quadrupole excitations 
are below 1 MeV, i.e. about a factor of 10 less than the 
energy of the dipole excitations. Since the quadrupole 
motion is ten times slower than the dipole one we use the 
adiabatic approximation: By means of RPA, we calcu- 
late the dipole absorption cross section <j 7 (E, f3 n , j n ) for 
a set of instantaneous deformation parameters {f3 n ,~f n } 
of the mean field. We find the probability P(/3„,7„) of 
each shape being present in the ground state by means 
of IBA-1, and obtain the total cross section as the inco- 
herent sum of the instantaneous ones, 



(1) 



In other words, we assume that the quadrupole deforma- 
tion does not change during the excitation of the nucleus 
by the absorbed photons, which sample the instantaneous 
shapes of the nucleus in the ground state. Accordingly 
we suggest the acronym ISS-RPA for the approach. 
Ref. ll| studied a collective dipole vibration coupled 



adiabatic approximation rather well reproduces the ex- 
act spreading. The ISS expression JT]) is of very general 
nature. All versions of RPA based on a deformed mean 
field could be used to calculate cr 7 (E, (3 n , 7„). Likewise, 
any model describing the the collective quadrupole mode 
could be used to obtain P(/3 n ,7„). 

We adopt^the quasiparticle version of RPA (QRPA) de- 
scribed in [lit ], which combines a triaxial potential with 
separable interactions. We changed to a Woods-Saxon 
potential with "Universal Parameters" . The pairing gaps 
are adjusted to the even-odd mass differences. The El- 
part of <r 7 is calculated with an isovector dipole-dipole 
interaction, the strength of which is adjusted to the ex- 
perimental position of the maximum of the GDR. The 
Ml-part is calculated with a repulsive isovector spin-spin 
and an isoscalar quadrupole-quadrupole interaction (for 
details cf. [l3| ). The cr~,{E) are calculated using the 
strength function method with a resolution of 100 keV. 

We describe collective quadrupole mode by means of 
IBA-1 [l4|, which is known to well reproduce the devel- 
opment of energies and E2-transition probabilities from 
spherical to well deformed nuclei through the transition 
region. We use the "extended consistent Q formalism" . 
The Hamiltonian and the E2-transition operators are 
given in fl5j ] . The parameters £, x, eg of the model are 
listed in Table HI For practical reasons we fix the boson 
number to Nb — 10, which turned out to be a sufficient 
flexible basis. The probability distribution P(/3„, y n ) is 
generated by means of the method suggested in [16[ . We 
consider the two scalar operators 



q 2 = [Q x <8 Q x 
qs = [Q x <8> [Q x <8> 
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which are formed by angular momentum coupling from 
the IB A quadrupole operators Q*. The two commuting 
operators q 2 and q^ are diagonalized in the space of states 
generated by coupling 10 bosons to zero angular momen- 
tum. The resulting eigenvalues q2, n and q3, n and the 
eigenstates |n) are linked to the deformation parameters 
as follows 
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to a harmonic quadrupole vibration and showed that the 



(4) 

The probabilities P(J3 n ,j n ) = \(0f\ n)\ 2 are the projec- 
tions of the IB A ground state |Q^~ ) on the set \n). The 
method is easily generalized to the case of shape coexis- 
tence. In this case one has two families of states, a and 
b. They are described by two sets of IBA parameters and 
the mixing coefficients c\ and c\ — 1 — c 2 -, which give the 
fraction of each family in the ground state. We generate 
P(flna,lna) and P{f3 nbl "/ nb ) separately. The total distri- 
bution is then given by c 2 P{/3 n v,lnv), u = a,b, and the 
sum in ([1]) runs over all combinations {n, v}. 

We applied the method to the nuclides 88 Sr, 90 Zr and 
92 ~ 100 Mo, for which the combination of earlier (7, n) 
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FIG. 2: (Color online) Photo-absorption cross section of 
94 Mo. MM: RPA for the equilibrium deformation, ISS: 
RPA averaged over the probability distributions for shapes in 
Fig. HI ISS-CD: ISS folded with a Lorentzian of width aE 2 . 
Experimental data from [13] 



measurements [ll with recent (7,7') experiments at the 



ELBE facility [B, 0, El HI provided absorption cross 



sections <J-y{E) covering the whole energy range from 
the GDR down to few MeV. The IBA parameters were 
obtained by fitting the energies and B(E2) values of 
the lowest 0+, 2+, 4+ states taken from the ENSDF 
data base and from 2l|, |22| . We set the boson number 
Nb = 10, but otherwise followed Refs. [2l|,[22j|, where we 
took the mixing compositions from. Fig. [T] demonstrates 
that the resulting instantaneous shapes are widely dis- 
tributed across the (3 — 7 plane, which reflects the tran- 
sitional nature of the considered nuclei. In the cases of 
90 Zr and 98,100 Mo the two coexisting families of shapes 
are clearly recognized (blue and red). 

Figs. [2j OH and |4] compare the results with the experi- 
mental data. We include the RPA results for the equilib- 
rium deformations (labelled by MM in the figures) listed 
in Tab. [H which were calculated by means of the micro- 
macro method [H, Let us first consider 4fMo 52 
shown in Figs. [U and [3l It has a spherical equilibrium 
shape, but pronounced transitional character, which is 
demonstrated by the wide distribution of shapes in Fig.[T] 
The spherical RPA (MM) shows the expected strong fluc- 
tuations of ct 7 , reflecting the degeneracy of the spherical 
single particle levels. The substantial Landau fragmen- 
tation generates dipole strength in the threshold region 
which is of the order of the experimental one. Sampling 
the different instantaneous shapes (ISS) largely wipes out 
the strong spikes, because the spherical degeneracy is 
lifted. It shifts additional strength into the region below 
14 MeV, such that tr 7 comes close to the data. At the cen- 
ter of the GDR, the ISS peak is broader and lower than 
the spherical RPA one. However, it is still too narrow 
and too high because the collisional damping is missing 
Q. We include it by folding the ISS cross section with 
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FIG. 3: As Fig. [2] for the energy range around the neutron 
emission threshold, where S n =9.68, 9.15, 8.64.8.29 MeV for 
N = 52, 54, 56, 58, respectively. Data from 



a Lorentzian of width Y = aE 2 , which depends on the 
photon energy as expected for collisional damping 0, [|| . 
The coefficient a is chosen to reproduce the experimental 
07 at the maximum of the GDR, which gives a =0.0105 
MeV- 1 for the neutron number N =50 and 0.014 MeV" 1 
for N > 50. The resulting curve labelled by ISS-CD re- 
produces the data very well. The agreement of ISS-CD 
with the data is as good as in Fig. [5] for the other nuclides 
not shown on this scale. 

Figs. [3] and 2] zoom into the region near the nucleon 
emission thresholds. ISS well describes the experimental 
<7 7 on the average, while still fluctuating too much. The 
fluctations are wiped out in the ISS-CD curve, which 
reproduces the data very well. The inclusion of colli- 
sional damping (ISS-CD vs. ISS) hardly shifts additional 
strength from the GDR to the considered energy region. 
Hence, the dipole strength near the nucleon emission 
threshold is determined by the Landau fragmentation of 
the instantaneous shapes, each of which contributes with 
its probability being present in the ground state. 

In order to quantify the contributions to the width 
of the GDR peak, we folded the MM and ISS results 
for 4IM052 with a Lorentzian of constant width T, cho- 
sen such that the peak height agreed with the exper- 
iment. For the MM case we found T=4.6MeV and 
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FIG. 4: As Fig. [3] The mean of the proton and neutron 
thresholds are (S p + S n )/2 =10.86, 10.16, 10.06 MeV for Z = 
38, 40, 42 , respectively. Data from [TJ, dj 



tended in a systematic way or it can be replaced by any 
large-amplitude description for the collective shape de- 
grees of freedom. 

In summary, we propose a novel method (ISS-RPA) 
for calculating the dipole strength function of nuclei with 
large-amplitude shape fluctuations, which combines the 
Interacting Boson Model (IBA) with the Random Phase 
Approximation (RPA). The method bases on the exis- 
tence of two time scales: the slow shape dynamics and 
the fast dipole vibrations. Instantaneous Shape Sampling 
(ISS) assumes that the photo-absorption occurs at a fixed 
shape, described by RPA, the probability of which given 
by IBA. We studied the nuclides with (Z=38, 40, 42, 
^=50) and (Z=42, N=52, 54, 56, 58). ISS-RPA very 
well reproduces the experimental photo-absorption cross 
sections <7 7 . Around the nucleon emission thresholds, <r 7 
is determined by the Landau fragmentation and the fluc- 
tuating shapes. ISS-RPA may be used for calculating 
<7 7 of unstable nuclei needed in astrophysical simulations 
of violent stellar events for a firm understanding of the 
origin of chemical elements in the universe. 

This work was supported by the DFG project 
KA2519/1-1 (Germany) and the US DOE grant DE- 
FG02-95ER4093. 



for ISS r=4.0 MeV. Since the experimental value is 
MeV, the contribution of the Landau fragmen- 
tation is r LF =(5. 7-4.6) MeV= LlMeV. ISS contributes 
r 7SS =(5. 7-1. 1-4.0) MeV=0.6 MeV. The strongest con- 
tribution Tcd— 4.0 MeV results from collisional damp- 
ing. Hence for the "spherical" nucleus 4IM052 the dis- 
tribution of dipole strength near the peak of the GDR is 
dominated by the collisional damping whereas Landau 
fragmentation and shape fluctuations dominate in the 
low-energy tail of the GDR. The contribution of shape 
fluctations to the peak width increases with N, going over 
to the splitting into two peaks characterizing a well de- 
formed shape. Applying the same procedure to the other 
Mo-isotopes, we found for the collisional damping width 
rcD=(4.0, 3.6, 3.5, 3.9) MeV and for the combined con- 
tribution from the fluctuating shape and Landau frag- 
mentation r LF+IS s=0- 7 i 2 - 7 > 2 - 5 > 4 -°) MeV for ^=94- 
100, respectively. 

We reproduce the experimental cr 7 near the nucleon 
emission thresholds using input parameters that are fixed 
by independent information, namely the location of the 
GDR and low-lying quadrupole excitations. This holds 
the promise that the concept of ISS-RPA will improve the 
prediction of the dipole strength for the unstable nuclides 
passed in the stellar events. We chose the combination, 
IBA-1 and RPA with dipole-dipole interaction, because 
it is simple and the modules were available. Our sim- 
ple RPA can be replaced by more sophisticated versions 
based on modern density functionals, which are expected 
to provide reliable dipole strength functions for unstable 
nuclei. Likewise, the IBA-1 phenomenology can be ex- 
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